Effects of dilution and disorder on magnetism in diluted spin systems 
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The influence of configurational disorder on the magnetic properties of diluted Heisenberg spin 
systems is studied with regard to the ferromagnetic stability of diluted magnetic semiconductors. 
The equation of motion of the magnon Green's function is decoupled by Tyablikov approximation. 
With supercell approach, the concentrations of magnetic ions are determined by the size of the 
supercell in which there is only one magnetic ion per supercell in our method. In order to distinguish 
the influence of dilution and disorder, there are two kinds of supercells being used: the diluted and 
ordered case and the diluted and disordered case. The configurational averaging of magnon Green 
function due to disorder is treated in the augmented space formalism. The random exchange integrals 
between two supercells are treated as a matrix. The obtained magnon spectral densities are used 
to calculate the temperature dependence of magnetization and Curie temperature. The results are 
shown as following: (i) dilution leads to increasing the averaged distance of two magnetic ions, further 
decreases the effective exchange integrals and is main reason to reduce Curie temperature; (ii) spatial 
position disorder of magnetic ions results in the dispersions of the exchange integrals between two 
supercells and slightly changes ferromagnetic transition temperature; (iii) the exponential damping 
of distance dependence obviously reduces Curie temperature and should be set carefully in any 
phenomenological model. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.50.Pp, 85.75.-d 
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I. INTRODUCTION 

Following the discovery of a ferromagnetic transition at 
temperatures in excess of 100K, the diluted III-V mag- 
netic semiconductors (DMS), which are realized by dop- 
ing a semiconducting host material with low concentra- 
tions of magnetic impurities (typically manganese), have 
attracted a great deal of attention from both the experi- 
mental and theoretical point of view due to their poten- 
tial in spintronics applicationsi*2iMiSiSi2i2iS. 

In DMS, low concentrations of magnetic impurities car- 
rying localized magnetic moments (spins) form a diluted 
spin system. The random spatial distribution of the mag- 
netic impurities breaks the translational symmetry of the 
crystal and thus greatly complicates the theoretical de- 
scription of the materialiSiiiiiSii 3 .. Several different theo- 
retical methods 14 ' 15,16 ' 1 - have been performed to get the 
transition temperatures. In some stages of these calcu- 
lations, disorder effects have been completely neglected 
or treated within the mean-field approximation (MFA). 
Monte Carlo (MC) simulationa 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 seem 
to provide a better way to include the positional 
disorder, but these theories usually assume classical 
spins. A proper treatment of the positional disorder 
of the localized moments and their quantum nature is 
necessarj*i2i2&. Recently, ab initio calculations with su- 
percell approachesSZiS&SSiSi are used to investigate the 
effects of disorder. Alternatively, the stochastic series ex- 
pansion (SSE) quantum MC (QMC) method with L x L 
supercell is used to investigate the order-disorder tran- 
sition in the diluted two-dimensional Heisenberg model 
with random site dilution^. But there are still doubts 
concerning the effect of disorder on magnetism. 



In fact, the effect of disorder on magnetism is an 
old and important problem in diluted spin system o 32 ! 33 , 
although only during the last few years have disorder 
effects in DMS been considered. In these early pa- 
pers, a typical method is using coherent potential ap- 
proximation (CPA) that is initially developed by Soven 
and TayloD 3 ^ 5 . to treat the dynamics of a random 
Heisenberg-type Hamiltonian. However, most of these 
methods only include the short-range interaction and are 
hard to extend to the long-range exchange interaction. 
Up to now, lots of methods of investigating spin sys- 
tems such as SSE-QMC 3 i and density matrix renormal- 
ization group thoery 3 ^ only take into account short-range 
or nearest-neighbor interactions. But the experimental 
fact of the transition temperature with low concentra- 
tions implies that the exchange interaction between the 
localized moments is long-ranged in DMSi* 3 . It means 
the long-range exchange interactions should be also in- 
cluded in the theoretical calculation of the magnetization 
and the transition temperature. 

Another method for calculating the properties of 
disordered systems is the augmented space formalism 
(ASF) 3 ^ 3 ^, which is introduced by Mookerjee 3 ^ and cen- 
ters around averaging functions of independent random 
variables. Rather than expanding the Green function 
in some manner and then averaging an appropriate set 
of terms, the random problem is transformed into an 
ordered one which is defined in a larger Hilbert space. 
This new Hilbert space is referred to as the augmented 
space, which may be described as the direct product of 
the Hilbert space spanned by the original Hamiltonian 
with a " disorder" space that describes the random vari- 
ables. On transforming to this augmented space, a new 
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nonrandom Hamiltonian can be defined such that config- 
urational averages in real space for the random solid are 
equal to inner products in the augmented spac o 4 ?' 41 ' 42 . 

In this article, we present a new approach, which com- 
bines the supercell approach and the ASF, to study the- 
oretically the influence of position disorder of magnetic 
ions on magnetization and the transition temperature in 
diluted spin systems on a disordered Heisenberg model. 
Firstly, the size of the supercell determines the concentra- 
tion of magnetic ions in the host materials. There is one 
magnetic ion per supercell which, however, can only oc- 
cupy the same site in a supercell (the diluted and ordered 
case) or can occupy any site in a supercell (the diluted and 
disordered case). In the diluted and disordered case, the 
distance between two magnetic ions therefore becomes 
a random variable and then the effective Heisenberg ex- 
change integrals, which are assumed to be a function of 
distance only, are random variables. In the framework 
of the ASF, the random exchange integrals are extended 
to matrices. Furthermore, the obtained spectral densi- 
ties are used to calculate the temperature dependence of 
magnetization and Curie temperature of systems. Signifi- 
cantly, the spins are treated quantum mechanically in our 
approach although we use the supercell approach just like 
MC simulations. Moreover, the long-range exchange in- 
tegrals are included. It should be mentioned here that the 
direct numerical diagonalization of the Green function 12 
is only applicable to the case of the finite size systems 
although it is a good method to include the long-range 
exchange integrals and to treat the spins quantum me- 
chanically. 

The article is organized as follows. The theoretical 
methods are described in section [H] Section Hill is con- 
cerned with the numerical studies and discussion. In sec- 
tion we conclude the article with a summary. 
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FIG. 1: A two dimension view of the sc {2 x 2 x 1} supercell 
(x — 25%) for (a) diluted and ordered case, and (b) diluted 
and disordered case. The small crosses represent the lattice 
sites and the black circles show the lattice sites occupied by 
magnetic ions. 



systems, we use the effective Heisenberg Hamiltonian 

N N 

H = -J2 J v S» • Si - ^gj^B^St • (1) 

i,j'=l i=l 

Here i and j label the sites occupied by the magnetic ions 
and Sj = (Sf, Sf , Sf) is the magnetic moment at lattice 
site i with lattice vector r^. Moreover, in our supercell 
approximation, % and j can be also referred to the labels 
of the supercell because there is only one magnetic ion in 
each cell. The exchange parameter is assumed to be a 
function of the distance |r^ — | only, independent of local 
environment. The Hamiltonian also contains a Zeeman 
coupling with external magnetic field B = (0, 0, B). 
Introducing the retarded magnon Green's function 



G ij {E) = {{St;Sj)y* 
where Sf = Sf ± iSf , its equation of motion reads 
(E - gj^B) Gij(E) = 2h 2 8 r] (Sf) - 2hJ2^m 



(2) 



II. THE MODEL 

For a diluted spin material Ai^B^ (A: non-magnetic; 
B: magnetic), the concentration x (x £ (0,1]) of mag- 
netic ion B is the ratio of the number of magnetic ions to 
all ions. If we consider a supercell that is built by non- 
magnetic ions with only one substituted by magnetic 
ions, the size of the supercell decides the concentration 
x. For example, for simple cubic (sc) systems, the con- 
centration x equals to 1/(1 x mx n) for the supercell with 
size {I x m x n}, where I, m and n refer to the numbers 
of ions in X, Y and Z axis, respectively. In Fig. as an 
example, we give two kinds of possible realizations of a 
sc {2 x 2 x 1} supercell that corresponds to the concen- 
tration x = 25%: (a) diluted and ordered case, and (b) 
diluted and disordered case. Each magnetic ion is located 
in the same lattice site of the supercell in (a), but in any 
possible site of the supercell in (b). 

To study the magnetic properties of diluted magnetic 



X Sj )) T E ~ ((SmSiiSj )) B ) ■ (3) 

To decouple the higher-order Green's function, we make 
the Tyablikov approximation on the rhs. of @. After 
rearranging, the equation of motion for Gij (E) is 



E 
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(4) 



where w — (E — g j [i-qB) / (2h (S z )) . Thus, the matrix 
form of Green's function can be expressed as: 

-l 



G(u) = h(u)I-K)~ 



(5) 



where I is the identity matrix and the matrix H has the 
elements Hij = (5y X^=i Jin — -hj which belong to the 
Hilbert space "K. 

In the concentrated case, after Fourier transforming, 
one can evaluate the Green's function from its momen- 
tum representation: 

G(<L " ) = ,-[J(fl)-J(q)l ■ (6) 
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Here J(q) = Yli Jri z Hi(*\) where J^i corresponds to 
the summation over the J-th shell with a distance Ri 
from a given site 0, zi is the total number of sites in the 
J-th shell and 7/(q) = j^J2r r e_ir "' q where the sum- 
mation over runs over each site located in the 7-th 
shell. 

The equation Jf)J) can also be applied to the calcula- 
tion of the diluted and ordered case if one considers all 
the distances r^- between two magnetic ions as equal to 
the Rjj between the respective supercells. Due to the 
translational symmetry of the supercell, it is chosen as 
the Wigner-Seitz unit cell. The shell structure is deter- 
mined by the spatial distribution of the supercells. For 
example, in Fig. ^ if the supercell i is set as the central 
supercell that refers to the 0-th shell, the supercells j and 
k belong to the first shell and the supercell I belongs to 
the second shell. 

But for the diluted spin system, equation (JHJ is not 
appropriate to be used directly because each lattice has 
different neighboring environment and the exchange inte- 
gral between two supercells is different. The translational 
symmetry on lattice and supercell is absent. According to 
(b) , the diluted and disordered case of the Fig. ^ the dis- 
tance between two magnetic ions, which can occupy any 
site in their own supercell, is a random variable, e.g. ry. 
Then all possible distances between two magnetic ions in 
different supercells are described as a set of independent 
random variables. For example, in the framework of the 
shell structure, the set of all possible distances between 
two magnetic ions, which belong to the n-th supercell of 
the 7-th shell and the central supercell, is {r%, r%, ... , r^} 
with the probability {c ri , c r2 , c rk } where c ri = 1, 
respectively. It leads to the exchange integrals JR In rep- 
resented as a random variables set { J ri , J T2 , . .., J Tk } with 
the probability density 



k 

i=i 



Jn) 



(7) 



It should be mentioned that, according to the shell struc- 
ture, the set of exchange integrals Jr t „ and the probabil- 
ity density p k {J Rln ) in the 7-th shell are the same when 
n runs over each supercell in the 7-th shell. So, one sim- 
plifies them as Jr s and p k (jR I ). 

For each random variable Jr : , according to the aug- 
mented space theorem, one can introduce a new Hilbert 
space (j) k such that p k (jRj) corresponds to a suitably cho- 
sen operator M k on 4> k of rank k, spanned by k possible 
values of Jr t . If |/q ) (usually 1,0, • • • ,0)) belongs to an 
orthonormal basis in (f> k , then 

p k {J Rl ) = -- i Hm Imif^Jh-M')- 1 ^) (8) 

7T ,/-» J Rl +lQ+ 

where Ik is the identity operator on the space 4> k . In other 
words, |/q ) and M k are chosen such that the spectral 
density of the operator M k with respect to |/q) is the 
given probability distribution. With the suitable choice 



of the basis, the tridiagonal matrix representation of M k 
has ai (I — 1, • • • , k) down the diagonal and b m (m = 
1, ■ • ■ , k — 1) down the off-diagonal positions, which can 
be obtained by the continued-fraction expansion: 
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p k (J Rl ) = — 



7T J- 



lim Im 



(9) 
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Thus the random exchange integral Jr x is transformed 
into the matrix representation M k in the Hilbert space 
4> k where the randomness of the exchange parameters is 
completely included. 

Now, in the framework of the ASF, each supercell has 
the same spin term and the exchange integral from the 
central supercell to one special shell has the same matrix 
form. No supercell in the configurational distribution is 
distinguishable from any other. In other words, in the 
diluted spin system, the translational invariance that is 
destroyed on lattice is recovered on the supercell while 
the randomness is masked in the exchange integral ma- 
trix M k . Note, the exchange integral matrix M k is the 
same within one and the same shell, but different among 
the different shells. Here, we want to mention an inter- 
esting technique: the combinatorial method^ that re- 
duces initial eigenproblem to a less-dimensional one. In 
the case of finite systems, considering symmetry prop- 
erties including methods of algebraic combinatorics, the 
Hamiltonian matrix is transformed to a block (quasidiag- 
onal) form and the numerical solutions of eigenproblem 
are high precision and small resultant. In our approaches, 
due to recovering in the translational invariances, one can 
use the equation JHJ to study an infinite system. 

In the ASF, the configurational averages in real space 
are replaced by inner products in the augmented space 
E = 'K ® where the "disorder" space $ = (j) 1 ® (f> 2 <g> 
■ • • <8 (j) 1 ® • • • with (p 1 = Y[V <E>4> kl being zj times the 
direct product of 4> kl and zi is the total number of the 
supercells in the 7-th shell. Moreover, \Fq) = \Fq) 
IT^o ) (g> • • • <g> \Fq ) (g> • • • is an orthonormal basis in the 
"disorder" space $, where |F 7 ) = Iff ®|/o x )- Note, in 
the expression of If^ 1 }, we add the index 7 to represent 
|/o) in the /-shell. About the details of constructing 
the augmented space for a randomly disordered system 
with independent site-occupation variables, we refer the 
reader to a series of papers2242iii4^. The elements of 
Green function can be expressed as 

GonH = (7? ®7^o|?i(^I-h) 1 17?„ ® F ) (10a) 

) • (10b) 

ij m 

Here, if i refers to the central supercell and j refers to 
the n-th supercell that is located in the 7-th shell, then 



J« = I 1 



r/-l 



(gi M 7 <g> I 



i+i 



, (11a) 



M 1 = (H®I kl )®M k ' ®(]J®I kl ) , (lib) 



71+1 
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where Yim> ®Ifej means (m — m' + 1) times the direct 
product of l kl and I L = fj* 1, l kL is the direct product of 
all zl identity matrices of rank in the L-th shell. 

Considering the translational symmetry of the super- 
cell, after Fourier transforming, one can evaluate the av- 
eraged Green's function: 



G 00 (q,w) = (F \h[wl-H{q) 
H(q) = Q(0) - Q(q) , 

where, in the shell structure form, 

Q(q) = £Q'(q) , 



\Fo) 



Q 7 (q) 

M 7 (q) 



i 1 



rJ-1 



iM J (q)(g)I 



1+1 



(12a) 
(12b) 

(13a) 

,(13b) 
(13c) 



EMi(q), 

n=l 

n—1 zj 



n+l 



where Mr 1 (q) = Mj 7/(q) (the proof being given in 
Appendix) . 

Furthermore, the averaged Green function (G(u>)) can 
be evaluated from its momentum representation 



(GH)=^G 00 (q, 



(14) 



where the Lambin-Vigneron algorithm 4 ^ is used to do 
summation in q-space over the Brillouin zone. 

Then, the magnon spectral function can be expressed 

as 

A(u) = --lm(GH) . (15) 
Using the Callen equation, magnetization reads 

. (1 + S + *) * 2S+1 + (S - (1 + *) 2S+1 



{S*) = fi- 



(1 + * 



,2S+1 



$2S+1 



(16) 

where the average magnon number 'J can be calculated 
by 



dw- 



,2h{s*)ui/k B T _ i 



(17) 



Considering (S z ) — > in the limit T — ► Tc, one can 
get the transition temperature 



knT, 



B±C 



\hS{S + l) ^ 



(18) 



III. NUMERICAL STUDIES 

In this section, the influence of dilution and disorder 
on magnetization and Curie temperature is investigated 
in diluted spin systems. For simplicity, we only consider 
the case of a three-dimensional system on a simple cubic 
lattice with the lattice constant a — 1 and an infinitesi- 
mal external magnetic field. In addition, the shrink of the 
q-space is included because of the Fourier transforming 
based on the supercells. 

To study the influence of the range of a ferromag- 
netic exchange interaction on magnetization, respec- 
tively, power-law, RKKY and damped-RKKY type ex- 
change interactions are used: 



' j /(— ) 4 

./„l/?) = < Jo M?) - (f)cos(f )] /(f ) 4 , 

y e-w[ sin (f)-(f)cos(f)]/(£)4 



(19) 



where a is the lattice constant, Jo is the effective nearest- 
neighbor interaction strength and R is the distance be- 
tween two magnetic ions. If the damping factor Rq — > oo, 
the damped-RKKY exchange interaction is equal to the 
RKKY exchange interaction. In order to calculate prac- 
ticably, one needs to do cutting about the exchange inter- 
action at a distance R cut - ff- Considering the Fourier 
transformation based on the supercells and comparing 
the calculations of the diluted and ordered case and of 
the diluted and disordered case, the diluted and ordered 
case is used to decide the shell structure of systems. If the 
distance between one supercell and the central supercell 
is larger than R cut - ff, the exchange interactions will 
be considered as zero. It should be mentioned that the 
equation l|19|l is used just for illustration of the present 
approach, although the neglected effect of disorder on the 
value of exchange integrals 23 ' 24 ^ should be included. In 
addition, superexchange effects and the influence of vir- 
tual bound states on the value of exchange integrals^^ 
are also neglected. 

Fig. shows the magnon spectral density of the di- 
luted and ordered case and the diluted and disordered case 
for several concentrations x that correspond to the dif- 
ferent sizes of the supercells. Because of the cut-off dis- 
tance R cu t-of f = 2.0a, the exchange parameters between 
the nearest-neighbor supercells are only included for the 
{2x2x2} supercells' structure that corresponds to the 
concentration x — 12.5%. That is the reason why the 
calculation of the {2x2x2} structure in the diluted and 
ordered case shows the well-known symmetric shape of a 
simple cubic density of states except for the change of the 
bandwidth. For the {1x1x2} and {2x2x1} supercells' 
structures, in the diluted and ordered case, there are more 
than one peak because of including the exchange inter- 
actions belonging to next nearest-neighbor (nnn) or other 
longer-range supercells. In addition, there is one peak in 
the {1x1x2} and {2x2x1} structures being in the same 
position as that in the {2x2x2} structure. It means all 
the magnetic ions and the exchange interactions in the 
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FIG. 2: The magnon spectral density for R cut ^ f f/a= 2.0 on 
an sc lattice for various concentrations x (that correspond to 
the different sizes of the supercells) for the diluted and ordered 
case and the diluted and disordered case with the effective 
nearest-neighbor integral Jo = 1.0, S = 2.5 and the damping 
factor Ro = 2.0a. 



{2x2x2} structure are also included in the {1x1x2} and 
{2 x 2 x 1} structures. It should be noted here that the 
long-range exchange interactions between magnetic ions 
are included in the nearest-neighbor supercells because 
the exchange interactions are based on the positions of 
magnetic ions. Fig. El also shows that the peaks of the 
magnon spectral density move toward lower energies for 
decreasing concentrations of the magnetic ions. It shows 
that dilution increases the magnon spectral density for 
lower energies at the cost of the magnon spectral density 
at higher energies. 

Disorder also enhances the low energy part of the 
magnon spectral density, especially for the {1x1x2} 
and {2 x 2 x 1} structures, if one compares the diluted 
and ordered case and the diluted and disordered case in 
Fig. |21 However, it is clearly shown in the figures that 
the influence of dilution is more important than that of 
disorder. Furthermore, the obviously separated peaks in 
the diluted and ordered case disappear in the diluted and 
disordered case because the random position distribution 
of the magnetic ions leads to a dispersion of the exchange 
integrals, while the exchange integral between two super- 
cells is a single value in the diluted and ordered case. At 
the same time, it is clear that the dispersion of the ex- 
change integrals leads to the long tails that extend to 
high energy magnons' area, especially for the {2x2x2} 
structure in the panel of power-law interactions. 

It is quite interesting to note that the damping fac- 
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FIG. 3: Magnetization (S z ) as function of temperature T for 
Rcut-of f /a — 2.0 on an sc lattice for various concentrations 
for the diluted and ordered case and the diluted and disordered 
case with Jo = 1.0, S = 2.5 and Ro = 2.0a. 

tor also increases the magnon spectral density for lower 
energies. From the expression of the damped-RKKY ex- 
change integrals, one finds that the exponential damping 
will decrease the strength of the short-range interaction 
less than that of the long-range one. Compared with 
the RKKY exchange interaction, it leads to the clearer 
separation of the exchange integral distribution for the 
damped-RKKY exchange interaction in the diluted and 
ordered case. It is especially clear for the {2x2x1} 
structure. But in the diluted and disordered case, the 
changing of the separation is not so obvious due to the 
dispersions of the exchange integrals. 

The temperature dependence of magnetization, which 
is calculated from the magnon spectral density in Fig. |21 
is shown in Fig. |3 The figures show that all curves 
of magnetization decreases to zero monotonically and 
smoothly with increasing temperature from T = to 
T = T c . From the equation (|18fl. one knows that increas- 
ing magnon spectral density at lower energies will lead to 
decrease the Curie temperature. It is clearly shown by 
the influence of dilution on Curie temperature because 
dilution strongly shifts the magnon spectral density to 
lower energies. Dilution increases the distance between 
two magnetic ions and further decreases the strength of 
the exchange integrals. In other words, the influence of 
dilution on Curie temperature mainly comes from de- 
creasing of the effective exchange interactions. 

Disorder also influences Curie temperature T c , but not 
so important as dilution. The random distribution of 
the magnetic ions leads to the dispersions of the ex- 
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change integrals. That increases the lower energy part 
of magnon spectral density and is the reason why disor- 
der lowers the Curie temperature. It is especially clear 
for the {1x1x2} and {2x2x1} structures. However, for 
the {2x2x2} structure of power-law interactions, Curie 
temperature T c is slightly higher in the disordered case 
than that in the ordered case because the dispersions also 
increases the high energy part of magnon spectral den- 
sity. The final effect of disorder on T c is the results of 
competition between the high energy magnons and the 
low energy magnons. If one uses an average value or 
a special value to substitute the dispersion of the ex- 
change integrals, just like the mean-field approximation 
(k,BT c ~ J2q J e ^( < l)) or the diluted and ordered case 
in our calculation, the influence of the dispersions of the 
exchange integrals will be washed out. 

The results show that the Curie temperature is mainly 
influenced by dilution or the effective exchange interac- 
tions. It may explain the following fact. For a long time, 
it was rather surprising that the models that neglect the 
disorder and treat the exchange interaction in MFA pro- 
vide a Curie temperature in a good agreement with the 
experiment. Although the positions of magnetic ions is 
randomly distributed, one may get a reasonably calcu- 
lated Curie temperature, which is near to the experimen- 
tal Curie temperature, by adjusting suitably the effective 
exchange integrals. Moreover, given the level of dilution 
of the magnetic ions, the results suggest that the way to 
get high Curie temperature DMS materials is increasing 
the effective exchange integrals Jo- 

Comparing the calculation of RKKY and damped- 
RKKY, we find there is a big influence of the exponential 
damping on Curie temperature, which is even a more im- 
portant factor than the position disorder of magnetic ions 
for Rq = 2.0a in our calculation. The reasons to intro- 
duce the exponential damping is due to substitutional 
disorder or to the half-metallic character of (Ga,Mn)As 
and (Ge,Mn) alloy, eto^i^S^. In our theory, only the 
latter is meant because the positional disorder is already 
explicitly taken into account in our calculation. The obvi- 
ous influence of damping on Curie temperature suggests 
that the damping factor should be chosen carefully if one 
only wants to incorporate the damping exponential item 
to include the positional disorder influence, especially in 
model calculation. 



IV. SUMMARY 

The aim of this article is to study the influence of di- 
lution and disorder on the ferromagnetic properties of 
diluted spin systems. By combining the supercell ap- 
proach and the ASF, we study theoretically the influence 
of disorder on magnetization and Curie temperature in 
diluted spin systems on a disordered Heisenberg model. 
Firstly, the size of the supercell is used to determine the 
concentration of magnetic ions because there is only one 
magnetic ion per supercell in our calculations. In order 



to investigate the influence of dilution and disorder on 
Curie temperature, the positions of magnetic ions in the 
supercells can be arranged in two ways: (i) the diluted 
and ordered case in which the magnetic ion is only able 
to occupy the same lattice within each supercell and (ii) 
the diluted and disordered case in which the magnetic ion 
is able to occupy any lattice point in the supercell. For 
the diluted and ordered case and the concentrated case, 
the equations of motion of the magnon Green's function, 
which is decoupled by making Tyablikov approximation, 
can be solved directly from Fourier transformation be- 
cause of the translational symmetry of the supercells. For 
the diluted and disordered case, the effective Heisenberg 
exchange integrals, which are assumed to be a function of 
distance only, are random variables because the positions 
of two magnetic ions in the supercells is randomly dis- 
tributed. By using ASF, the random exchange integrals 
between two supercells are extended to matrices, there- 
with restalling translational symmetry, in order to cal- 
culate the averaged magnon spectral density. Then, the 
obtained averaged spectral densities have been used to 
calculate the temperature dependence of magnetization 
and the Curie temperature. Significantly, the long-range 
exchange integrals are included and all the calculations 
correspond to an infinite system and there is no cutting 
in the real space. 

The resulting theory is then solved numerically for a 
simple cubic system. The long-range ferromagnetic ex- 
change integrals including power-law decaying, RKKY 
type and damped-RKKY type are used to calculate 
the temperature dependence of magnetization and Curie 
temperature. The results show that dilution and the 
damping exponential item increase the magnon spectral 
density for lower energies at the cost of the magnon spec- 
tral density at higher energies. That leads to a decrease 
of the Curie temperature of spin systems. The influence 
of dilution on Curie temperature mainly comes from de- 
creasing the effective exchange interactions. The effect 
of disorder on Curie temperature is more complicated. 
The random distribution of the magnetic ions leads to 
dispersions of the exchange integrals, which mainly de- 
crease Curie temperature but sometimes increase T c in 
the concentration area of our calculations. The role of the 
exponential damping of distance is obvious and should 
be set carefully in any phenomenological model. From 
our calculations, to attack the problem of high temper- 
ature ferromagnetic DMS materials, the effective way is 
increasing the effective exchange integrals. 
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V. APPENDIX In addition, for 



{1)2, •••,zr} — > {R/2, • • ■ i R/21 1 R-ji} j 



For Ry, if the index i refers to the J-th shell and j 
refers to the j-th supercells in 7-th shell, we write Fourier 
transformation of Jy = J(i2y) as: : 

Q(q) = y2e^ iRir ' l J(R i j) {l: 2 ,---,z/} — > {R/ Z/ , R/i, • • • , R/ (Z7 _2) ; , 

fly 

_ [ e -' R n qj(i? n ) -| 1- e~ 4Rli i' q J(i?i Zl )] correspondingly, one can write 

+ [e-' iR21 - q J(i? 2 i) + • • • + e- lR ^3(R 2z2 )] 

+ ■■■ M J (q) = e -** R » (Mj 1 ® J fcl <g> 7 fc , ® • ■ • ® 7 fc , ) 

+ [e-^^'^J^n) + • • • + e-'^-iJ^)] + e -^- R «( /fe/ ® m*' ® J fcj ® • •• ® I fcl ) 

+ ••■ + ••• 

= J2 Q'W + e-** Rw (4, ® 4, 8) • • • ® 4, ® ) , 

where is the total number of the supercells in the 7-th 
shell and 



Q 7 (q) - J2 e lRln '^(Rm) M / (q) = e -*i-R^(Af* J <»J fcl (8/ fa ®-.-<8)J fcl ) 

n=l 



n=i To sum zj expressions of M 1 (q) and divided by zj, one 



If we refer the n-th supercells to that of the central po- 
sition R/„, i.e., 



can get 



{l,2,...,(z J -l),z ; }^{R /ll R ;2 ,..,R J , 2l _ 1) ,R j2; }, M J (q) = ^M^(q) 

M J (q) can be expressed as 



71=1 

n— l 



Ml(q) = (n®Ife I )®^(q)®(n®Ifcx), 
M J (q) = e- ,€ »- R "(M^(8)/ fc/ ®/ fc7 <8>-- ■«/*,) i „+i 



+ e -**«-« g> MY ® I fcj ® • • • ® 4,) 

+ ■■■ where Mj J (q) = M^7/(q) and 77 (q) 

+ e- iq ' R ^ (J fel ® J fc , ® • • • ® ® M*') i E^=i e-* R '"- q . 
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